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Abstract 



Hard photon emitted from energetic heavy ion cohisions is of very interesting since it does 
not experience the late-stage nuclear interaction, therefore it is useful to explore the early-stage 
information of matter phase. In this work, we have presented a first calculation of azimuthal 
asymmetry, characterized by directed transverse fiow parameter F and elliptic asymmetry co- 
efficient V2, for proton- neutron bremsstrahlung hard photons in intermediate energy heavy-ion 
collisions. The positive F and negative V2 of direct photons are illustrated and they seem to be 
anti-correlated to the corresponding free proton's flow. 
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The properties of nuclear matter at differ- 
ent temperatures or densities, especially 
the derivation of the Equation-of-State 
(EOS) of nuclear matter, are one of the 
foremost challenges of modern heavy-ion 
physics. Since heavy ion collisions provide 
up to now the unique means to form and 
investigate hot and dense nuclear matter 
in the laboratory, many experimental and 
theoretical efforts are under way towards 
this direction. Because of their relatively 
high emission rates, nucleons, mesons, light 
ions and intermediate mass fragments, 
produced and emitted in the reactions, are 
conveniently used to obtain information on 
the reaction dynamics of energetic heavy 
ion collisions. However, these probes in- 
teract strongly with the nuclear medium 
such that the information they convey may 
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bring a blurred image of their source. For- 
tunately, energetic photons offer an attrac- 
tive alternative to the hadronic probes [1 . 
Photons interacting only weakly through 
the electromagnetic force with the nuclear 
medium are not subjected to distortions 
by the final state (neither Coulomb nor 
strong) interactions. They therefore de- 
liver an undistorted picture of the emitting 
source. For hard photons, defined as 7- 
rays with energies above 30Mel^ in this 
paper, many experimental facts supported 
by model calculations [1,2,3] indicate that 
in intermediate energy heavy-ion colli- 
sions they are mainly emitted during the 
first instants of the reaction in incoherent 
proton-neutron bremsstrahlung collisions, 
p + n— i>p + ?T, + 7, occurring within the 
participant zone. This part of hard pho- 
tons are called as direct photon. Direct 
hard photons have thus been exploited to 
probe the pre-equilibrium conditions pre- 
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vailing in the initial high- density phase of 
the reaction [4,5]. Aside from the dominant 
production of hard photons in first-chance 
p-n collisions, a significant hard-photon 
production in a later stage of heavy-ion 
reactions, called as thermal photons, are 
also predicted by the Boltzmann-Uehling- 
Uhlenbeck (BUU) theory [6,7]. These ther- 
mal photons are emitted from a nearly 
thermalized source and still originate from 
bremsstrahlung production by individual 
p-n collisions, which was also confirmed by 
the experiments at last decade [8,9]. 

In this work, we take the BUU trans- 
port model improved by Bauer [10]. The 
isospin dependence was incorporated into 
the model through the initialization and 
the nuclear mean field. The nuclear mean 
field U including isospin symmetry terms 
is parameterized as 
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where po is the normal nuclear matter den- 
sity; p, pn, and Pp are the nucleon, neutron 
and proton densities, respectively; equals 
1 or -1 for neutrons and protons, respec- 
tively; The coefficients a, b and a are pa- 
rameters for nuclear equation of state. Csym 
is the symmetry energy strength due to the 
density difference of neutrons and protons 
in nuclear medium, which is important for 
asymmetry nuclear matter (here Cgym = 32 
MeV is used), but it is trivial for the sym- 
metric system studied in the present work. 

For the calculation of the elementary 
double-differential hard photon production 
cross sections on the basis of individual 
proton-neutron bremsstrahlung, the hard- 
sphere collision was adopted from Ref. [11], 
and modified as in Ref. [12] to allow for en- 
ergy conservation. The double differential 
probability is given by 
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Here R is the radius of the sphere, a is 
the fine structure constant. Pi and /3/ are 
the initial and final velocity of the proton 
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Fig. 1. (a) Time evolution of hard photon pro- 
duction rate for the reaction ^^Ca Ca col- 
lisions at 30 MeV /nucleon for semi-central 
events (40-60%). (b) Time evolution of re- 
duced maximum density Pmax/po (closed cir- 
cles) and reduced average density (p) / po (open 
circles) of the whole reaction system in the 
same reaction. The blue dashed line represents 
the time when the system ends up till the first 
expansion stage, and in the panel (a) it sep- 
arates direct photons (on the left side) and 
thermal photons (on the right side). 

in the proton-neutron center of mass sys- 
tem, and 9^ is the angle between incident 
proton direction and photon emitting di- 
rection. More details for the model can be 
found in Ref. [10]. 

In this paper, we simulate the reaction of 
^°Ca + Ca coUisions at 30M eV /nucleon, 
and use the £'05' with the compressibility 
K of 235 MeV {a = -218 MeV, b = 164 
MeV, a = 4/3) for the nuclear mean field 
U. As a first attempt to extract the pho- 
ton's azimuthal asymmetry, we only take 
the semi-central events (40 — 60%) as an 
example in this Letter. 

In Fig. 1 we show the time evolution of 
production rate of bremsstrahlung hard 
photons as well as the time evolution of 
system densities, including both maximum 
(closed circles) and average density (open 
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circles). We found that hard-photon pro- 
duction is sensitive to the density oscilla- 
tions of both the maximum and the average 
density during the whole reaction evolu- 
tion. When the density of collision system 
increases, that is in the compression stage, 
the system produces more hard photons. 
In contrary, when the system expends, the 
hard photon production decreases. Actu- 
ally, the density oscillations of the colliding 
heavy ions systems can be observed in the 
experiments via hard-photon interferome- 
try measurements [13,6]. Apparently, hard 
photons are mostly produced at the early 
stage of the reaction. Combining the time 
evolution of the nuclear density, we know 
that this part of hard photons are domi- 
nantly emitted from the stage of the first 
compression and expansion of the system. 
Thereafter we call these photons, emitted 
before the time of the first maximum ex- 
pansion of the system (t = 80 fm/c in this 
reaction), as direct photons (on the left 
side of blue dashed line in Fig. 1(a)). It is 
also coincident with the definition of direct 
photons above. And we call the residual 
hard photons produced in the later stage 
as thermal photons (on the right side of 
blue dashed line in Fig. 1(a)). So in the 
simulation, we can identify the produced 
photon as direct or thermal photon by the 
emitting time. Because of the sensitivity 
to the density oscillations of coUiding sys- 
tem, hard photon may be sensitive to the 
nuclear incompressibility [6,7]. 

It is well known that collective flow is 
an important observable in heavy ion 
collisions and it can bring some essen- 
tial information of the nuclear matter, 
such as the nuclear equation of state 
[14,15,16,17,18,19,20,21,22,23]. Anisotropic 
flow is defined as the different n — fh har- 
monic coefficient Vn of the Fourier expan- 
sion for the particle invariant azimuthal 
distribution [15]: 

dN °° 

— ocl + 2^v„cos(n(/)), (3) 

where is the azimuthal angle between the 
transverse momentum of the particle and 
the reaction plane. Note that the z-axis 



is defined as the direction along the beam 
and the impact parameter axis is labelled 
as X-axis. Anisotropic flows generally de- 
pend on both particle transverse momen- 
tum and rapidity, and for a given rapidity 
the anisotropic flows at transverse momen- 
tum pt {pt — ^pI +pI) can be evaluated 
according to 

Vn{Pt) = {cos{n(j))), (4) 

where (■ ■ ■ ) denotes average over the az- 
imuthal distribution of particles with trans- 
verse momentum pt, Px and Py are projec- 
tions of particle transverse momentum in 
and perpendicular to the reaction plane, 
respectively. The flrst harmonic coefficient 
Vi is called directed fiow parameter. The 
second harmonic coefficient f 2 is called the 
elliptic flow parameter V2, which measures 
the eccentricity of the particle distribution 
in the momentum space. 

In relativistic heavy-ion collisions az- 
imuthal asymmetry of hard photons have 
been recently reported in the experiments 
and theoretical calculations [24,25,26,27]. 
It shows a very useful tool to explore the 
properties of hot dense matter. However, 
so far there is still neither experimental 
data nor theoretical prediction on the az- 
imuthal asymmetry of hard photons in 
intermediate energy heavy ion collisions. 
Does the direct photon also exist azimuthal 
asymmetry so that it leads to non-zero di- 
rected transverse flow or elliptic asymme- 
try parameters in the intermediate energy 
range? Moreover we know that direct pho- 
tons mostly originate from bremsstrahlung 
produced in individual proton-neutron col- 
lisions, and free nucleons are also emitted 
from nucleon-nucleon collisions. Does the 
azimuthal asymmetry of the direct photons 
correlate with the one of free nucleons? To 
answer the above question, we focus on 
the azimuthal asymmetry analysis for both 
photons and protons in this Letter. 

Fig. 2 shows the time evolution of the di- 
rected flow parameter vi and elliptic flow 
parameter V2 for hard photons and free 
protons. Before we take further calculation 
and explanation, people should be cau- 
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tious about the word of "flow" for photons. 
Since flow is associated with collectivity 
caused by multiple interactions, which are 
exhibited by the nucleons, but not by the 
photons. The photon emission pattern is 
basically a result of the nucleon flow, and 
not a photon flow per se. However, How- 
ever, in order to compare the results be- 
tween photons and protons, we still called 
Vi as directed flow parameter and f 2 ellip- 
tic flow parameters for photons somewhere 
in texts. Considering the nearly symmet- 
ric behavior for directed flow parameter 
(fi) versus rapidity, here we calculate the 
average vi over only the positive rapidity 
range, which can be taken as a measure 
of the directed transverse flow parameter. 
For emitted protons (open circles in the 
figure) which are experimental measurable 
and are identified in our BUU calculation 
as those with local densities less than po/8, 
the onset of flows occurs around t = 30 
fm/c before that the system is mostly in 
fusion stage and protons are seldom emit- 
ted. The negative directed flow parameter 
Vi of free protons essentially stems from the 
attractive mean field. Up till t ~ 120 fm/c 
when the system is in the freeze-out stage, 
the directed flows become saturate. For 
the elliptic asymmetry parameter w 2 of free 
protons, the positive values indicate of the 
preferential in-plane emission driven by 
the rotational collective motion due to the 
attractive mean filed. Similarly, the elliptic 
asymmetry parameter becomes saturate 
in the freeze-out stage. However, there 
are obvious difference for proton-neutron 
bremsstrahlung photons (solid circles in 
the figure) in comparison to protons. Con- 
trary to the negative directed transverse 
flow and positive flow, directed photons 
shows the positive vi and the negative 
V2 before t = 80 fm/c, i.e. the azimuthal 
anisotropy is shifted by a phase of 7r/2. The 
times corresponding to the peak or valley 
values of flows roughly keep synchronized 
with the compression or expansion oscilla- 
tion of the system evolution. For the late- 
stage thermal photons after t = 170 fm/c 
the azimuthal asymmetry vanishes, i.e. vi 
and V2 fades-out. 

From the above calculations, we learn that 
thermal photons in the later stage of reac- 
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Fig. 2. The time evolution of V2 (a) and vi 
(b) for hard photons (closed circles) and free 
protons (open circles). 

tion are emitted from a more thermalized 
system, they prefer more isotropic emission 
(i.e. the vanishing "flow" parameters) than 
direct ones produced in the pre-equilibrium 
stage. Thereafter we only consider direct 
photons to discuss the azimuthal asymme- 
try results. For protons, we take the values 
of flows when the system has been already 
in the freeze-out time at 180 fm/c. 

The directed transverse flow parameter at 
mid-rapidity can be also defined by the 



slope: F = ^F"^' , where (y)cm 

is the rapidity of particles in the center of 
mass and {px) is the mean in-plane trans- 
verse momentum of photons or protons in 
a given rapidity region. In Fig. 3(a) and 
(b), we show {px) plotted versus the cm. 
rapidity Ucm. for direct photons (a) as well 
ci-s {px) plotted versus the reduced cm. 
rapidity {y/ybeam)c.m. for free protons (b). 
The errors shown are only statistical. A 
good linearity was seen in the mid-rapidity 
region (—0.5, 0.5) and the slope of a linear 
fit can be defined as the directed transverse 
flow parameter. The extracted value of the 
directed transverse flow of direct photons 
is about -1-3.7 MeV/c, and that of free pro- 
tons is about —12.4 MeV/c. Thus direct 
photons do exist the directed transverse 
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Fig. 3. (a) Average in-plane transverse mo- 
mentum of direct photons function of cm. 
rapidity for semi-central events (40 — 60%). 
The dashed line segment is a fit over the 
mid-rapidity region —0.5 < Ucm. < 0.5. (b) 
Same as the panel (a) but for free protons. The 
dashed line segment is a fit over the mid-ra- 
pidity region -0.5 < {y/ybeam)c.m. < 0.5. (c) 
and (d) are the azimuthal distributions of di- 
rect photons and free protons, respectively, 
and both of them are fitted to the 4i/i order 
Fourier expansion. 

asymmetry even though the absolute value 
is smaller than the proton's flow, and its 
sign is just opposite to that of free protons. 

As Eq. 3 shows, elliptic flow is defined as 
the second harmonic coefficient V2 of an 
azimuthal Fourier expansion of the par- 
ticle invariant distribution. In order to 
extract the value of elliptic asymmetry 
coefficient V2 and reduce the error of fits, 
we fit the azimuthal distribution to the 
Ath order Fourier expansion. Shown in 
Fig. 3(c) and (d), direct photons demon- 
strate out-of-plane enhancement and the 
V2 is about —2.7%. Whereas, for free pro- 
tons, azimuthal distribution displays the 
preferential in-plane emission and the f 2 is 
about +7.2%. Furthermore, we can extract 
the transverse momentum dependence of 
the elliptic asymmetry coefficient V2- Fig. 4 
shows V2 of direct photons (a) and free 
protons (b) function of transverse mo- 
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Fig. 4. V2 as a function of transverse momen- 
tum (pt) for direct photons (a) and free pro- 
tons (b). 

mentum px- Similar to the directed trans- 
verse flow parameter, the values of elliptic 
asymmetry coefficient f 2 of direct photons 
and free protons also have the opposite 
signs at this reaction energy, i.e. reflecting 
a different preferential transverse emission 
in the direction of out-of-plane or in-plane, 
respectively. Meanwhile, the absolute val- 
ues of V2 for photons are smaller than the 
proton's values as the behavior of trans- 
verse flow. Except the opposite sign, we see 
that both f 2 have similar tendency with 
the increase of pt, i.e., their absolute values 
increase at lower pt, and become gradually 
saturated, especially for direct photons. 

To explain the above anti-correlation of 
anisotropic emission between direct pho- 
tons and free protons, we should note that 
direct photons originate from the individ- 
ual proton-neutron collisions. As Eq. 2 
shows, we can roughly consider that in the 
individual proton-neutron center of mass 
system, in directions perpendicular to in- 
cident proton velocity, i.e. 9^ = 7i/2, the 
probability of hard photon production is 
much larger than that in the parallel di- 
rection, i.e. 6^ = 0, which is in agreement 
with the theoretical calculations and the 
experiments [29,30], that causes hard pho- 
ton preferential emission perpendicular to 
the motion plane of corresponding nucle- 
ons. As a whole, the azimuthal anisotropy 
of hard photons is shifted by a phase of 
7t/2 with respect to that associated with 
the anisotropy of nucleons, leading to the 
opposite signs of the values of F and V2 
between them. Consequently, azimuthal 
anisotropic emission of hard photon and 



5 



free nucleoli are anti-correlated, presenting 
the opposite behavior. 
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In conclusion, we have presented a first 
calculation of azimuthal asymmetry, both 
directed and elliptic asymmetry, for di- 
rect photons produced by proton-neutron 
bremsstrahlung from intermediate en- 
ergy heavy-ion collisions. It was, for the 
first time, presented that in the inter- 
mediate energy heavy-ion collisions the 
proton-neutron bremsstrahlung hard pho- 
ton shows non-zero directed transverse flow 
parameter and elliptic asymmetry coeffi- 
cient which have opposite sign to the corre- 
sponding free proton flow parameters. The 
time evolutions of azimuthal parameters Vi 
and V2 of hard photons exhibit rich struc- 
tures as the density oscillation of the sys- 
tem during the pre-equilibrium and ther- 
malization stage of reaction system. There- 
fore direct photons can server for a good 
probe to nuclear matter properties. Con- 
sidering that hard photons are dominantly 
produced by individual neutron-proton 
bremsstrahlung, so they are sensitive to the 
in-medium neutron-proton cross section, 
but not to the in-medium proton-proton or 
neutron-neutron cross section, that can be 
advantaged in the isospin dependent study 
of in-medium nucleon-nucleon cross sec- 
tion by direct photons. Of course, system- 
atic studies of the influences from equation 
of state, in-medium nucleon-nucleon cross 
section, impact parameter and incident 
energy etc on the azimuthal asymmetry of 
direct photon should be carried out. The 
progress along this line is underway. 
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